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Abstract. An analysis of the CKM matrix parameters within the /?fit approach is presented using 
updated input values with special emphasis on the recent sin2p measurements from BABAR and 
Belle. The QCD Factorisation Approach describing B kk,Kk decays has been implemented in 
the software package CKMfitter. Fits using branching ratios and CP asymmetries are discussed. 


STATISTICAL FRAMEWORK AND INPUTS 

In the Standard Model (SM) with three families, CP violation is generated by a single 
phase in the CKM matrix [p . This picture can be probed quantitatively by means of a 
global fit to all quantities sensitive to CKM elements in the SM. The analysis presented 
here is performed within the 7?fit statistical approach [Q], which is implemented in the 
software package CKMfitter [^. 

The quantity = —21ni: (yniod) L minimized in the fit, where the likelihood func¬ 
tion is defined by L{y ^od) = i^exp(jCexp --^theo(Tmod)) ■ i^theo(yQCD)- The experimental 
part, Lexp, depends on measurements, Xexp, and theoretical predictions, Xtheo, which are 
functions of model parameters, ymod- The theoretical part, i^theo, describes the knowl¬ 
edge on the QCD parameters, yqco ^ {ymod}, where the theoretical uncertainties are 
considered as allowed ranges. The agreement between data and the SM is gauged by 
the global minimum Xminymod’ determined by varying all model parameters ymod- For 
5^min vniod’ ^ confidence level (CL), expressing the goodness-of-fit, is computed by means 
of a Monte Carlo simulation. If the hypothesis “the CKM picture of the SM is correct” 
is accepted, CLs in parameter subspaces a, e.g. a = (p,fi) [0, are evaluated. For fixed 
a, one calculates A%^(a) = Xmin-fji^) ~ ^^minymod’ stands for all model parame¬ 

ters (including ygco) with the exception of a. The corresponding CL is obtained from 
CL(a) = Prob(Ax^(a), Adof), where Adof is the number of degrees of freedom, in general 
the dimension of the subspace a. Since the CL depends on the choice of the ranges for 
the ygcD, the results obtained in the fit have to be interpreted with care. 

The input values used in this analysis are listed in Tab. |I| For \Vub\, inclusive measure¬ 
ments from LEP and exclusive measurements from CLEO have been used. The prelimi¬ 
nary CEEO lepton endpoint analysis [Q] using moments obtained from B A^yis not yet 
included. Eor 114*1, inclusive measurements from EEP, the measurements of 5 —s- D*iv 
at zero-recoil and the moments analysis from CEEO [|^], using inclusive B —> Xc£\ and 
B —> XsJ decays, have been combined. The uncertainty on Ama has been significantly 
reduced due to the measurements from the 5-factories [j|]. However, the constraint on 
(p,fi) is not improved since it is dominated by the theoretical uncertainty on 



TABLE 1. Left; input values for the global fit. Right; ht results quoted for CL > 5 %. 


Input Parameter 

Value 

Fit Output 

Range (CL >5%) 

\Vud\ 

0.97394 ±0.00089 

X 

0.2221 ±0.0041 

\Vus\ 

0.2200 ±0.0025 

A 

0.763-0.905 

\Vcd\ 

0.224 ±0.014 

P 

0.07-0.37 

\VcA 

0.969 ±0.058 

fi 

0.26 - 0.49 

\Vub\ 

(3.49±0.24±0.55)-10-3 

7(10-3) 

2.2-3.7 

\Vcb\ 

(40.4 ±1.3 ±0.9)-10-3 

sin2a 

1 

o 

o 

1 

p 

Ui 

Anid 

(0.489 ± 0.008) ps-' 

sin2p 

0.59-0.88 

Anis 

Amplitude Spectrum 

Y 

37° - 80° 

\^k\ 

(2.271 ±0.017)-10-3 

Ams (ps-*) 

14.6-32.0 

sin2p 

0.793 ±0.102 

fB,V^d (MeV) 

192-284 

ot,(MS) 

(166 ±5) GeV 

Bk 

0.52-1.68 

OTe (GeV) 

1.30±0.10 

m, (GeV) 

95-405 

Bk 

0.87 ±0.06 ±0.13 

«(Az,^7t“vv)-(10-") 

1.6-4.4 

hoc 

L38±0.53 

«(A+ ^ 7t+vv) • (10-'^) 

4.8-9.4 

fB.V^d (GeV) 

(0.230 ± 0.028 ± 0.028) GeV 

«(B+ ^ T+Vt) • (10-3) 

5.6-23.8 

riB(MS) 

0.55 ±0.01 

1.16 ±0.03 ±0.05 

«(B+ ->P+V^) • (10-3) 

2.2-9.4 
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FIGURE 1. Left; constraints in the (p,r|) plane where p and r| are dehned by p = p(l — A-^/2), 
fj = T|(l — A-^/2). Shown are the 5 % CLs for the individual constraints as shaded areas and the lo- 
and 2c- contours from sin2p. In addition, the 5 % and 90 % CL contours for the combined ht are drawn. 
Right; constraints in the (p,f|) plane if the likelihood ratio is used for the Airig constraint. 


For Anis, the most recent combined amplitude spectrum from 0 is included in the fit us¬ 
ing a modified version of the standard amplitude method 0] . If the amplitude spectrum 
is translated into a likelihood ratio SB, a stronger constraint is obtained. However, to 
our knowledge, it has not been demonstrated so far that the likelihood ratio can be inter- 
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FIGURE 2. Left; Constraints in the (p,fi) plane using (see text) and the bound on R. The CL is 
indicated by the shaded areas; dark-grey: > 90 %, gray: 90 % — 32 %, light-gray: 32 % — 5 %. Right; 
Hypothetical constraints for summer 2002. 


preted as a probability density funetion. Henee, we use the more eonservative method of 
Ref. [0 for the numerical analysis presented here. For sin2p, the world average is used. 
It should be noted that the most precise measurements from BABAR and Belle [|^] differ 
presently by about two standard deviations. 


FIT RESULTS 

The global minimum of the CKM fit is found to be ^ = 2.3, resulting in a 

goodness-of-fit of 71%. It quantifies the excellent agreement between experimental data 
and the CKM picture of the SM. Fig. shows the (p,fi) plane. Drawn are 5% CL 
contours from the single constraints using Am^/, Am^, \Vub/Vcb\ and respectively, 
and the lo- and 2o-contours for the four-fold ambiguity on (3 from sin 2(3. Shown in 
addition are the contours for the combined fit including sin2(3. The statistical precision 
of the sin2(3 measurement already competes with the indirect, theoretically limited 
constraints. Fig. |T] (right) illustrates the improved constraint when using the likelihood 
ratio for Am^. Selected numerical results are summarized in Tab. [T|. 

CHARMLESS TWO BODY DECAYS 

Constraints on the angle y can be obtained from B decays. Based on color 

transparency arguments, theoretical calculations such as the QCD Factorisation Ap¬ 
proach (FA) dl^j and the QCD hard scattering approach [(111] have been developed. Re¬ 
cently, the FA has been implemented in CKMfitter. At present, it is premature to in- 



















fer reliable eonstraints on the basis of these ealeulations due to open theoretieal ques¬ 
tions [|^ [^, 1^ . Data from the 5-faetories are not yet preeise enough to probe the 
ealeulations in detail. Henee, all fit results are marked by an appropriate “R&D” logo. 
For this review, a global fit to 5 ^ nn.Kn branehing fraetions and direet CP asymme¬ 
tries measured in self-tagging B Kn deeays has been performed within the frame¬ 
work of the FA, where most reeent experimental results from BABAR [|r3|]. Belle [ |T^ ] 
and CLEO [[l5l] have been used. The numerous theoretical parameters are let free to vary 


within the ranges given in Ref. [|^. We find ^ = 2.0 and conclude that data are 

consistently described within the FA. The best FA fits are found at y 80° and are in 
agreement with the constraints from the standard fit. 

Using less theoretical assumptions, ratios of branching fractions can be formed to derive 
constraints in the (p,fi) plane. As an example, the CP-averaged ratio 


R = 


Xg± ® 
X50 


( 1 ) 


provides the bound R > sin^y which is independent of the strong phases [|T^. Unfortu¬ 
nately, the present world average R = I.OV+q}^ leads to weak constraints only owing to 
the tails of the experimental errors. The ratio 






measured to be R* = O.VO^Qjg, can be used to derive bounds in the (p,fi) plane [[TtI] . 
An important input for the theoretical prediction of R* is the tree-to-penguin ratio (P/T) 
which can be determined experimentally using the relation 


( 2 ) 


£3/2 =^th-tanec- — 

Jit 


fK /2 ■ « (R± ^ 7t±7tO) 




(3) 


where Rth stands for SU(3) breaking corrections estimated in the FA to be [ ]lli| ] Rth = 
0.98 ±0.05. The bound on R* can be translated into a prediction if additional information 
on the strong phases is inserted 10. Adopting the values for theoretical ranges quoted 
in Ref. [|^, one obtains the constraints shown in Fig. ^ At present, the constraints re¬ 
main rather weak due to the limited experimental precision. The slight deformation of 
the shape pattern around y ~ 90° is due to the bound on R. For summer 2002, an inte¬ 
grated luminosity of 100 fb ^ is expected to be collected by each experiment, BABAR 
and Belle. The experimental precision will then start to provide interesting constraints, 
as can be seen from Fig. obtained assuming the present central experimental values 
and appropriately rescaling their errors. However, the constraints would still rely on the 
validity of some theoretical assumptions not yet fully explored. 

Within the FA, the P/T ratio for B is predicted. Compared to the present exper¬ 

imental error on the time-dependent asymmetry Sim = O.OS^q j^ ±0.11 from BABAR, 
the quoted theoretical uncertainty is much smaller [j^]. In Fig. ^ (left), the constraints in 
(p,fi) from Sjirt are shown using P/T from FA where theoretical uncertainties have been 
neglected. The right plot shows the constraints when also using sin2(3. 
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FIGURE 3. Left; Constraints on (p,ri) from S-jm using the penguin-over-tree ratio from FA. Right; 
Constraints when using also sin2p. 
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